Limits of slow-light in photonic crystals 
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While ideal photonic crystals would support modes with a vanishing group velocity, state-of-the 
art structures have still only provided a slow-down by roughly two orders of magnitude. We find that 
the induced density of states caused by lifetime broadening of the electromagnetic modes results 
in the group velocity acquiring a finite value above zero at the band gap edges, while attaining 
superluminal values within the band gap. Simple scalings of the minimum and maximum group 
velocities with the imaginary part of the dielectric function or, equivalently, the linewidth of the 
broadened states, are presented. The results obtained are entirely general and may be applied to any 
effect which results in a broadening of the electromagnetic states, such as loss, disorder, finite-size 
effects, etc. This significantly limits the reduction in group velocity attainable via photonic crystals. 

PACS numbers: 42.70.Qs, 42.25.Bs, 78.67.-n 



The interest in slow-light phenomena has increased sig- 
nificantly in recent years, due to the potential applica- 
tions in areas such as optical processing [l|, 0, 0, 01 , quan- 



tum information processin, 
emission [7], and sensing 




, enhanced spontaneous 
Strongly dispersive peri- 



odic structures with dielectric functions that vary on the 
length scale of the wavelength of light can now be fab- 
ricated with impressive precision and resolution 1^, 11 1 . 
By exploiting the close analogies with condensed matter 
theory, these so-called photonic crystals can give rise to 
such novel ^tical phenomena as e.g., enhanced Raman 
scattering 12], increased stimulated emission ML super- 
prism behavior ^13.] . and negative refraction [l4|. Due 
to their highly dispersive nature, photonic crystals have 
emerged as excellent candidates for sources of slow-light 
phenomena [isl . [l"" 



17 



18l |. However, even though ideal 
structures would in principle support modes of vanishing 
group velocity, state-of-the art structures have still only 
provided a slow-down by roughly two orders of magni- 
tude [isl ITij l . The limits imposed on the minimum at- 
tainable group velocity in photonic crystals have been 
studied in various contexts, such as, e.g., fabrication dis- 
order J3, HO], lossy dielectrics [21], and finite-size ef- 
fects 1221. It is the aim of this paper to generalize these 
findings and to show that they may all be presented in 
the context of broadening of electromagnetic modes and 
the resulting induced density of states. 

The existence of photonic band gaps in photonic crys- 
tals is accompanied by van Hove singularities in the den- 
sity of states (DOS) near the band gap edge, which results 
in a vanishing group velocity. In ID the group velocity 
is inversely proportional to the density of states, while 
in higher-dimensional structures it is inversely propor- 
tional to the projected one-dimensional density of states 
along the propagation direction, as in mesoscopic elec- 



tron transport. A formal proof can be given in terms 
of Wigner-Smith group delay formalism, thus also ap- 
plying to inhomogcneous structures of finite extension 
(See e.g. Ref. 23] and references therein). Letting go 
denote the one-dimensional projected density of states 
along the propagation direction of a homogenous mate- 
rial with a dielectric function given as the average value of 
the dielectric function of the photonic crystal, we identify 
three different regimes of interest: (i) a long wavelength 
regime, where the properties of the photonic crystal are 
independent of the detailed geometrical composition, and 
thus g{Lu) ~ go, (ii) a slow-light regime, where g{u!) > go, 
and (Hi) a "superluminal" regime, where g{uj) < go- We 
stress that 17(0;) in all cases refers to the one-dimensional 
projected density of states along the propagation direc- 
tion, and that indeed this is the case whenever we speak 
of the density of states. Exp eriments have confirmed the 



ITI . |24| ] and superlumi- 



existence of both slow [iSl . [if 
nal [2^ regimes. 

Consider the dispersion relation near the edge of a 
photonic band gap. The exact shape of the dispersion 
relation naturally depends on the geometry and dimen- 
sionality of the photonic crystal. However, common to 
all band structures is a vanishing slope at the band gap 
edge. If we expand the band structure near the band gap 
edge uj(K)=LUo we thus obtain 



uj{k) ~ Wo + a{k - K)^ 



(1) 



where a is the group velocity dispersion (GVD). Here we 
have ignored any higher-order contributions to the ex- 
pansion. Higher-order dispersion is of vital importance 
for the propagation of pulses in the slow-light regime [13] , 
but we ignore it for now as we arc only interested in 
studying the actual value of the group velocity. The 
group velocity becomes 
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duj 



Re( ^) = Re(2y^ 



wo)) 



(2) 
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illustrating the square-root divergence of the density of 
states at the band gap edge. Let us consider the ef- 
fect of a finite imaginary part of the dielectric function 
e=e' + ie" of either of the constituents of the photonic 
crystal. To capture the effect of a small imaginary part 
we apply standard electromagnetic perturbation theory 
as described in detail in Ref. |27|. The first order term 
in a perturbative expansion in the imaginary part of the 
dielectric function is 



Alu 



(EUe"|E; 



Vi 



2 (E| e' |E) ' 



(3) 



where the integral in the denominator is restricted to 
the region Vi of the constituent containing the imaginary 
dielectric part. The effect is then an imaginary shift of 
the frequency 



where 



f 



(E|6ME)^^ 
(Ele' IE) 



(4) 



(5) 



is the fraction of dielectric energy localized in the corre- 
sponding dielectric [2^. Making the substitution luq 
ujQ + AujQ in Eq. ^ this means that at the band edge 
the group velocity becomes 



afuQ — 



(6) 



Despite its simplicity, this results is very general and in- 
dependent of the particular photonic crystal geometry, 
be it of one, two or three-dimensional nature. The result 
indicates that the vanishing group velocity and the corre- 
sponding divergence of the density of states is resolved in 
the case of a finite imaginary part of the dielectric func- 
tion. Furthermore, the group velocity has a sub-linear 
dependence on e", so the reduction in the group index is 
significant, even for small imaginary parts of the dielec- 
tric function. In deriving this result we have of course as- 
sumed that the slow-down is achieved at the band edge, 
where ideally Vg = 0. Other schemes exist for achieving 
slow-light within a flattened band, wherein the group ve- 
locity is non-zero [sij. However, including such a term in 
the previous derivation merely adds an additional term 
to Eq. (O inside the square root, and as such simply 
serves to increase the lower limit of the attainable group 
velocity. Due to the effects on pulse broadening, much 
attention has been devoted to finding structures with low 
group velocity dispersion in the slow-light regime 2^, 3^ . 
Interestingly, our analysis reveals that such structures 
have the additional benefit of reducing the minimum at- 
tainable group velocity. Indeed, in the limit of vanishing 
GVD our analysis shows that the group velocity attains 
the ideal value of zero, provided that any higher-order 



dispersion is negligible. Alternatively, we may consider 
the effect of e" as a finite broadening of the order e"uj of 
the electromagnetic modes Thus, any effect caus- 

ing such lifetime broadening of the modes is subject to a 
similar analysis as just described, and will thus increase 
the minimum group velocity as the square-root of the 
finite linewidth of the broadened modes. This very gen- 
eral analysis applies to any photonic crystal geometry, 
and any effect of lifetime broadening, such as e.g., loss, 
disorder, or finite-size effects. Of course, the analysis is 
based on a perturbative approach and as such is only 
applicable in the case of small perturbations of a per- 
fect photonic crystal. In the limit of strong disorder, for 
example, Anderson location will significantly alter the 
delay time statistics of the photonic crystal S^, 11]. In 
this analysis we limit ourselves to the case where the 
mean-free path is sufficiently long so that the concept of 
group velocity is meaningful. In passing, we note that 
another effect of the broadening of the modes is that a 
density of states is induced within the band gap. The 
induced density of states mid-gap is ^pbg oc e"/(Aijj)^, 
where Au is the width of the gap 27| . The result of this 
is that superluminal group velocities are attained within 
the band gap, wherein the perturbative analysis reveals 
that Vg oc 1/e". 

In the following we support these general findings by 
two illustrative examples, where the source of the broad- 
ening is a finite conductivity of the dielectric solid of 
an air-dielectric structure. We first consider the simple 
case of a Bragg stack with layers of width ai=0.8A and 
02=0. 2A for air and dielectric, respectively, where A is 
the lattice constant. We take ei=l for air, and assume 
that the dielectric solid can be described by a frequency 
independent real part £2=9 of the relative dielectric func- 
tion. The imaginary part is modeled as e2=cr / {sqi^) , 
where a is the conductivity of the dielectric and eo is 
the vacuum permittivity. This is similar to the Drude 
model for the optical response of metals, and is chosen 
so that the Kramers-Kronig relations are fulfilled. 

In Fig. [T] we show the band structure and correspond- 
ing density of states of the Bragg stack for the two cases 
of zero and nonzero conductivity of the dielectric solid, 
calculated using exact analytical expressions 3J| . At zero 
conductivity we note that two photonic band gaps exist, 
where there are no states with finite real /c-values. These 
PBGs are caused by the complete destructive interfer- 
ence of transmitted and reflected waves at each layer of 
the Bragg stack. The loss caused by the nonzero conduc- 
tivity makes the destructive interference incomplete, and 
the PBGs are no longer fully developed. Consequently, 
a density of states is induced in the band gap, and the 
divergence of the density of states at the band gap edge 
is resolved. The effect of this is two-fold. Firstly, the 
group velocity at the band edge acquires a finite value 
larger than zero because now g^oo, and secondly, the 
group velocity inside the PBG attains superluminal val- 
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FIG. 1: (a) Photonic band structure of the Bragg stack. The 
real parts of the fe-vectors are shown in blue, while the imag- 
inary parts are shown in red. The dashed lines indicate the 
case where the conductivity of the dielectric solid is nonzero. 
The photonic band gaps are indicated by yellow shading, (b), 
The corresponding density of states, with the case of nonzero 
conductivity shown in red. The dashed horizontal line in- 
dicates the density of states of a homogenous material with 
e = (aiei + a2e'2)/A. 
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ues larger than c because in this region g<go. We stress 
that while the group velocity in the superluminal regime 
may exceed the speed of light in vacuum, this implies no 
loss of causality, but is rather due to a strong reshaping 
of the optical pulse, as discussed for example by Biittiker 
and Washburn 35] . 

In Fig. [2^ we show the group velocity as a function 
of the normalized frequency, for three different values of 
the dimensionless parameter (;=aA/{2Trceo), character- 
izing the loss of the dielectric solid. To calculate the 
group velocities we use the definition Vg = duj/d(Re[k]), 
which yields the same results as the formal definition 
Vg = Re{duj / dk) [11]. Even for a low value of <j=0.01, we 
find that the group velocity no longer drops below ap- 
proximately 1 /80th of the speed of light in vacuum, and 
for ^=0.1 the minimum group velocity is one tenth of the 
speed of light in vacuum. Because of the lower overlap of 
the electromagnetic field with the dielectric solid in the 
air bands, the effect of loss on the group velocity is less 
pronounced for these bands. While the value of the min- 
imum group velocity is clearly highly dependent on the 
conductivity, the frequency ujq at which the group veloc- 
ity is at its minimum remains nearly constant, except for 
a very slight red-shift at high values of the conductivity. 
In Fig. [2)3 we show the maximum and minimum values 
of the group velocity as a function of the dimensionless 
conductivity. The maximum values of the group veloc- 
ity follow very precisely the l/e2 dependence expected 



FIG. 2: (a) Group velocity as a function of frequency 
for three different values of the dimensionless parameter 
C=(TA/(27rceo), characterizing the conductivity of the dielec- 
tric solid. The solid black line corresponds to the case of zero 
conductivity. The dashed horizontal line indicates the group 
velocity of a homogenous material with e — (aiei -I- 02 £2) /A. 
(b) The minimum group velocities of the dielectric and air 
bands of the first PBG, and the maximum group velocities of 
the first and second PBG. All are shown as a function of the 
dimensionless conductivity. The dashed black lines indicate 
fits to the numerical data of the form A^/q for the minimum 
group velocities, and of the form A/q for the maximum group 
velocities. The inset shows the geometry of the Bragg stack. 



from our initial analysis. Assuming a similar distribu- 
tion of the induced DOS in both band gaps, this analysis 
also indicates why the maximum group velocity falls off 
more rapidly in the case of the second, narrower band 
gap. The minimum group velocity is found to follow al- 
most exactly the expected -^/ej dependence. We note 
again that the dependence of the minimum group veloc- 
ity on loss is less pronounced for the air band, due to 
the smaller overlap of the electromagnetic mode with the 
lossy dielectric. We have verified that these simple scal- 
ing laws of the minimum and maximum group velocities 
hold for any parameters of the Bragg stack, and also in 
the case of oblique incidence. 

We next consider a two-dimensional photonic crystal 
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FIG. 3: (a) Group velocity as a function of frequency 
for three different values of the dimensionless parameter 
?=(TA/(27rceo), characterizing the conductivity of the dielec- 
tric solid. The solid black line corresponds to the case of 
zero conductivity. The dashed horizontal line indicates the 
group velocity of a homogenous material with e = ei + {e'2 — 
ei)7rr'^/A^. (b) The minimum and maximum group velocities 
as a function of the dimensionless conductivity. The dashed 
black lines indicate fits to the numerical data of the form 
for the minimum group velocities, and of the form A/q for the 
maximum group velocities. The inset shows the geometry of 
the photonic crystal. 



consisting of dielectric cylinders of diameter d/A — 0.55 
surrounded by air and arranged in a square lattice with 
lattice constant A. We take e'2 = (3.88)^ for the dielectric 
cylinders, corresponding to silicon, and model the loss via 
the Drude model, as for the Bragg stack. We consider 
normal incidence of TM modes. As for the Bragg stack, 
the effect of loss is to induce a density of states in the 
PBGs of the structure. The dispersion relations are cal- 
culated using the transfer matrix method [sj ■ In Fig. [3^ 
we show the group velocity as a function of normalized 
frequency for increasing values of the dimensionless con- 
ductivity. The results shown are for frequencies near the 
first TM photonic band gap. Similar trends are seen for 
the 2D photonic crystal as for the Bragg stack, namely 
that the density of states induced by the finite conductiv- 



ity of the dielectric resolves the divergence of the density 
of states at the band gap edge, thus causing Ug > at this 
point. Also, superluminal group velocities are attained 
inside the photonic band gap. In Fig. [31d the minimum 
and maximum group velocities are shown as functions of 
the dimensionless conductivity. As for the Bragg stack, 
the maximum and minimum group velocities are approx- 
imated very precisely by l/e2 and i/e^ dependencies, re- 
spectively, supporting our conclusion that these scalings 
are indeed universal for any geometry, and even dimen- 
sionality, of the photonic crystal. 

In this paper, we have for simplicity focused on con- 
ductivity as the cause of the smearing of the density of 
states. However, it is clear that any effect leading to a 
finite broadening of the electromagnetic states will result 
in a similar scaling of the minimum attainable group ve- 
locity with the order of the broadening. This includes 
loss, disorder, finite-size effects, etc. and consequently 
imposes significant limits on the minimum group veloc- 
ities attainable via photonic crystals. Even in the ab- 
sence of other broadening mechanisms, absorption will 
inevitably be present and thus represents the ultimate 
limiting factor on the attainable reduction in group ve- 
locity. In closing, we note that while we have focused our 
attention to photonic crystals, the analysis presented in 
this paper applies equally well to any slow-light scheme 
based on band structure effects, such as for example cou- 
pled resonator waveguides, where it has been shown that 
lattice disorder may severely limit the attainable reduc- 
tion in group velocity [s^. 
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